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Abstract. The Macaulay2 package CharacteristicClasses provides commands 
for the computation of the topological Euler characteristic, the degrees of the 
Chern classes and the degrees of the Segre classes of a closed subscheme of 
complex projective space. The computations can be done both symbolically 
and numerically, the latter using an interface to Bertini. We provide some 
background of the implementation and show how to use the package with the 
help of examples. 



1. Introduction 

The Macaulay2 [5] package CharacteristicClasses computes degrees of Chern and 
Segre classes of complex projective schemes. It also provides a command computing 
the topological Euler characteristic. Recall that the Chern classes of a smooth 
scheme arc defined to be the Chern classes of the tangent bundle. Moreover, the 
Segre classes of a scheme X embedded in P™ are defined to be the Segre classes of 
the normal cone CxP" of the embedding. Note that the degree of a d-dimensional 
cycle class, i.e., the weighted sum of classes of d-dimensional subvarieties, is defined 
to be the weighted sum of the degrees of the subvarieties. 

The computations done by the package CharacteristicClasses are based on a 
number of related algorithms which reduce the problem to the computation of the 
degrees of residuals, i.e., components complementary to a certain scheme in the 
intersection of hypersurfaces containing that scheme. The algorithm computing 
Chern classes is described in [5], the one computing Segre classes for possibly sin- 
gular schemes in [7], and the computation of the topological Euler characteristic 
is described in [11] . The residuals can be computed either symbolically or numer- 
ically, which yields a symbolical and a numerical version of the algorithms. Both 
versions are implemented in the package, the numerical version via an interface to 
Bertini [3]. 

There are other ways of computing the topological Euler characteristic and de- 
grees of characteristic classes. The Macaulay2 command euler computes the topo- 
logical Euler characteristic of smooth projective varieties by computing the Hodge 
numbers. In [2J, Aluffi describes algorithms for the computation of degrees of char- 
acteristic classes and the topological Euler characteristic. These algorithms are 
implemented in Aluffi 's package CSM, which is not part of Macaulay2 but can be 
obtained at www. math, f su. edu/^aluf f i/CSM/CSM. html. The different implemen- 
tations complement each other, as is shown in [7J. 

This article describes version 0.2 of the package CharacteristicClasses, which 
can be obtained at www. math. su. se/^jost/CharacteristicClasses .m2. Earlier 
versions provide computations of degrees of Chern and Segre classes, but not of the 
topological Euler characteristic. 
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2. Computing degrees of Chern and Segre classes 

We describe how to make Chern and Segre classes of projective schemes com- 
putationally tractable by computing their degrees. The standard reference for the 
notions used in this section, Chern classes, Segre classes and Chow groups, is |5]. 

Let X be a fc-dimensional closed subscheme of P™, embedded by i: X P", 
and denote its Chow group by A* (X) = Ad(X). Characteristic classes of 

X are elements of the graded group A*(X), whose generators are generally hard 
to compute. Hence algorithms for the computation of characteristic classes focus 
on computing a coarser invariant, the degrees of the classes. Let a = J2i a i[^i] 
be a cycle class in Ad(X), given as the weighted sum of classes of d-dimensional 
subvarieties Vi- Then the degree deg(a) of a is defined to be the weighted sum 
J2i a i deg(t^) of the degrees of the varieties Vt, seen as subvarieties of P™. One can 
also consider the pushforward of the cycle a to the Chow group of P n , which is 
known to be A*(P") = Z[H] / (H n+1 ) , where H is the class of a general hyperplane. 
The two viewpoints are equivalent because i*{a) — deg{a)H n ~ d . Both points of 
view, degrees of the classes and pushforward to the Chow ring of P n , are supported 
by the package CharacteristicClasses. 

The Chern classes ci(X), . . . , Ck(X) of a smooth scheme X are by definition the 
Chern classes ci(Tx), ■ ■ ■ , Cfc(Tx) of the tangent bundle of X. The total Chern class 
is the sum 1 + ci(X) + . . . + Ck(X). Furthermore, the Segre classes si(X,¥ n ), . . . , 
Sk(X, P") of a possibly singular X are the Segre classes si(CxP n ), ■ ■ • , Sk(CxV n ) 
of the normal cone C^P™ of X in P™. If the embedding is regular, the normal 
cone is a vector bundle, called the normal bundle. The total Segre class is the sum 
l + «i(X,P n ) + ... + s fc (X,P n ). 

The package CharacteristicClasses provides the commands chernClass and segre- 
Class which compute the degrees of the Chern and Segre classes of a given closed 
subscheme of P" or, equivalently, the pushforward of the total Chern and Segre 
class to the Chow group of P™. The closed subscheme X is given by a number 
of generators of a homogeneous ideal in a polynomial ring. The main idea of the 
algorithm is to relate the degrees of the Chern and Segre classes to the degrees of 
so-called residuals. One chooses randomly a number of hypersurfaces containing 
the scheme X, where the number of hypersurfaces is at least the codimension of 
X. According to a Bertini-type theorem, they intersect in X and some possibly 
empty component of expected codimension, the residual. Its degree can be com- 
puted either symbolically or numerically. Symbolically one computes the saturation 
using Grobner basis techniques. The numerical computations can be done using 
software for the numerical solution of polynomial equation systems. The package 
CharacteristicClasses implements both the symbolic version of the algorithms and 
the numerical, using an interface to Bertini [3] for the latter. More details on the 
algorithms can be found in [5] and [7]. 

3. Computing the topological Euler characteristic 

The topological Euler characteristic of a complex projective variety is the Euler 
characteristic of the underlying topological space with the usual Euclidean topology, 
i.e., the alternating sum of the Betti numbers. For the computations with the 
package CharacteristicClasses, we use that the topological Euler characteristic is 
equal to the degree of the top Chern-Schwartz-MacPherson class. More generally, 
the package can also compute the degrees of all the Chern-Schwartz-MacPherson 
classes. It uses the algorithm described in [11 , which reduces the computation 
of the degrees of Chern-Schwartz-MacPherson classes to the computation of Segre 
classes of singular subvarieties. For a more detailed introduction to Chern-Schwartz- 
MacPherson classes we refer to [5] and [TTj . 
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4. Using the package CharacteristicClasses 

We present three examples which demonstrate how to use the package Char- 
acteristicClasses. Observe that this article describes version 0.2 of the package 
CharacteristicClasses, but only version 0.1 is included in a standard installation of 
version 1.5 of Macaulay2, the most recent version at the moment of writing this 
article. More detailed information on the commands provided by this package can 
be found in its documentation, which is shown by e.g. using the command viewHelp 
CharacteristicClasses. 

Example 1. We compute the degrees of the Chern classes of a toy example, the 
twisted cubic. We start by loading the package and defining the ideal of the twisted 
cubic, which is generated by the 2-by-2 minors of a 2-by-3 matrix. Then we use 
the command chernClassList to obtain a list of the degrees of the Chern classes 
of the twisted cubic. The twisted cubic Ct w has dimension 1, hence it only has 
one Chern class Ci(C tw ). The output of the command chernClassList is the list 
{deg(C tw ),deg( Cl (C t w))}. 
Macaulay2, version 1.5 

with packages: ConwayPolynomials , Elimination, IntegralClosure , 
LLLBases, PrimaryDecomposition, ReesAlgebra, 
TangentCone 

11 : loadPackage "CharacteristicClasses"; 

— loading configuration for package "CharacteristicClasses" from 
file /. . . / .Macaulay2/init-CharacteristicClasses .m2 

12 : R = QQ[x,y,z,w] ; 

13 : twistedCubic = minors (2, matrix{{x,y,z},{y,z,w}}) 

2 2 
o3 = ideal (- y + x*z, - y*z + x*w, - z + y*w) 

03 : Ideal of R 

14 : chernClassList twistedCubic 

04 = {3, 2} 
o4 : List 

The computations above confirm that the twisted cubic has degree 3. As the 
twisted cubic is smooth, the degree of the top Chern class equals its Euler character- 
istic x(Ctw), which is related to the genus g(Ct w ) by \ = 2—2g. As deg(ci(Ct w )) = 2 
by the computations above, they confirm that the genus is and the twisted cubic 
is a rational curve. 

As said in Section 2, computing the degrees of the Chern classes is equivalent to 
computing the push- forward of the total Chern class to the Chow ring of P™. The 
total Chern class of the twisted cubic is 1 + ci(C tw ), its pushforward to the Chow 
ring of projective space is deg(Ct w )# 2 +deg(ci(Ct w ))ff 3 , where H is the hyperplane 
class. The pushforward can be computed using the command chernClass. 
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15 : chernClass twistedCubic 

3 2 
o5 = 2H + 3H 

ZZ[H] 

05 : 

4 

H 

All computations can also be done numerically using an interface to Bertini. 
Version 1.3 of Bertini needs to be installed and the package must be configured 
correctly. For more information on the configuration, use viewHelp "configuring 
Bertini". Then by using the value Bertini for the option ResidualStrategy, the 
computations are done numerically using Bertini instead of Grobner basis compu- 
tations. 

16 : chernClassList (twistedCubic , ResidualStrategy=>Bertini) 

06 = {3, 2} 

06 : List 

Example 2. We proceed with an example for the computation of Segre classes. 
The Whitney umbrella is a singular surface in P 3 . Over the reals its looks like an 
umbrella and the singular locus like the handle of the umbrella. We compute the 
degree of the first Segre class s±(S, P") of the singular locus S. Observe that in 
the following the singular locus is represented in Macaulay2 not as an ideal, but 
as a projective variety. All commands in this package work for both ideals and 
projective varieties. 

±7 : whitney = ideal (x~2*w - y~2*z) ; 

07 : Ideal of R 

18 : handle = Proj singularLocus whitney; 

19 : segreClassList handle 
o9 = {1, 0} 

o9 : List 

We see that the singular locus is a curve with degree 1, hence a line, and the 
degree of its first Segre class is deg(si(5, P™)) = 0. According to [IJ, the degree of 
the first Segre class is deg(si(5, P")) = v — 2, where v is the number of pinch points 
of the surface. It follows that the Whitney umbrella has two pinch points. Only 
one of them is visible in the real patch {w ^ 0} in which the Whitney umbrella 
usually is drawn, it is the tip of the umbrella. 

Example 3. We continue with an example from algebraic statistics which uses com- 
putations of the topological Euler characteristic, Example 2.2.2 from [BJ. The theory 
behind the computations is described in more detail in [IT]. The random censor- 
ing model with two events is a statistical model implicitely described by the ideal 
(2p pip 2 + p\p2 + P\p\ - P%P\2 + P1P2P12) in the polynomial ring C\p ,Pi,P2,Pi2], 
where po, pi, P2 and p\2 describe the probabilities of two events to occur be- 
fore or after a third event. Given experimental data one would like to compute 
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the values for the probabilities po, Pi, P2 and P12 which describe the data best, 
by maximizing the likelihood function. However, the function may have several 
stationary points and methods such as the Newton method may only find a lo- 
cal maximum. Hence it makes sense to define the maximum likelihood degree as 
the number of critical points of the likelihood function, as was done in [3]. By a 
theorem of Huh [10] , for a large class of examples the maximum likelihood degree 
equals the signed topological Euler characteristic of a certain open subvariety of the 
model. In this case, it is the topological Euler characteristic of the open subvariety 
V(2p p 1 p 2 +p\p2 +P1P2-P0P12+P1P2P12) \ V(p pip 2 pi2(po +Pi +P2 +P12)), where 
P0P1P2P12 7^ means that no probability should be zero, and po +P1+P2 +P12 7^ 
means that the probabilities should sum up to 1. We compute the topological Euler 
characteristic using the command eulerChar together with the inclusion-exclusion 
principle. 

110 : S = QQ[pO,pl,p2,pl2] ; 

111 : randomCensoring = ideal (2*p0*pl*p2 + pl~2*p2 + pl*p2~2 - 

p0~2*pl2 + pl*p2*pl2) ; 

011 : Ideal of S 

112 : boundary = ideal ( p0*pl*p2*pl2* (p0+pl+p2+pl2) ) + 

randomCensor ing ; 

012 : Ideal of S 

113 : eulerChar randomCensoring 

013 = 5 

114 : eulerChar boundary 

014 = 2 

It follows that the topological Euler characteristic of V(2popip2 +P1P2 +P1P2 ~ 
P0P12 +P1P2P12) \ V(poPiP2Pi2(po +Pi +P2 +P12)) is 5 - 2 = 3. Hence the maximum 
likelihood degree of the random censoring model is 3, which confirms the result in 
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